Van der Waals coefficients for systems with ultracold polar alkali-metal molecules 
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A systematic study of the leading isotropic van der Waals coefficients for the alkali-metal atom + 
molecule and molecule + molecule systems is presented. Dipole moments and static and dynamic 
dipole polarizabilities are calculated employing high-level quantum chemistry calculations. The 
dispersion, induction, and rotational parts of the isotropic van der Waals coefficient are evaluated. 
The known van der Waals coefficients are then used to derive characteristics essential for simple 
models of the collisions involving the corresponding ultracold polar molecules. 
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I. INTRODUCTION 

Ultracold chemistry in sub-microkelvin regime has 
emerged as one of the most exciting fields in atomic and 
molecular physics [1-3]. By tuning magnetic field across 
Feshbach resonances one can combine two free atoms into 
a bound state, with binding energy of order of MHz and 
then, with elaborated laser techniques, coherently trans- 
fer them into the deeply bound states - including the ab- 
solute rotational- vibrational-electronic ground state. At 
present two alkali-metal dimers have been produced in 
this manner KRb and Cs2 [4, 5]. It is also worth men- 
tioning that LiCs molecules in the vibrational ground 
state have been produced by photoassociation followed 
by spontaneous emission [6]. At present, many experi- 
mental groups have focused on production of other hete- 
ronuclear alkali-metal dimers hoping to obtain ultracold 
quantum gases of polar molecules, stable with respect to 
the atom exchange and trimer formation [7] . Such quan- 
tum gases of polar molecules will be used to explore new 
ideas in quantum information theory [8, 9], quantum sim- 
ulations of condensed-phase physics [10], or fundamental 
studies of chemical reactions [11]. 

Description of chemical processes in sub-microkelvin 
regime is extremely difficult, because the full quan- 
tum calculation for such systems is nearly impossible. 
Thus only few quantum dynamics studies of ultracold 
atom-|-diatom collisions employing global potential en- 
ergy surfaces have been performed so far concentrating on 
the homonuclear spin-polarized systems [12-18], where 
single-electronic-state approach provides good approxi- 
mation [19]. Even then, the quantum dynamics calcula- 
tions for heavier system are very challenging and have not 
been yet performed despite the fact that the correspond- 
ing quartet potential energy surfaces are rather simple 
[20, 21]. 
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Theoretical treatment of non-spin-polarized systems 
would be even more challenging. The calculations of tri- 
atomic and tetraatomic interaction potentials in such a 
case would have to include many active electrons and cou- 
pled potential energy surfaces, which at present is very 
far from routine. The following quantum dynamics cal- 
culations, especially in the presence of external fields, 
would be extremely demanding. The interaction poten- 
tials involving alkali-metal atoms and dimers are likely 
to be strongly anisotropic, and therefore the basis sets 
for such calculations would have to be very large. On 
the other hand, there is a very small number of observ- 
ables as outcome of ultracold collisions. After all, in a 
laboratory we do not record state-resolved cross sections 
but only loss rates from the state prepared before the 
experiment. Thus, the recent theories of ultracold colli- 
sions [22-25], formulated to explain current experiments 
in this field, use only few simple parameters that catch 
the essential physics. Importantly enough, the feature of 
the intermolecular interaction that matters the most is 
the long-range shape of the interaction potential, usually 
represented analytically by the well-known van der Waals 
expansion with the most important term —CqR^^ (R is 
the distance between the monomer centres of mass). 

Properties of the alkali-metal dimers have been inten- 
sively studied using electronic structure methods. A sys- 
tematic study of the dipole moments of all possible alkali- 
metal dimers was published by Aymar and Dulieu [26], 
and Deiglmeyr et al. [27] reported a systematic study of 
the static dipole polarizabilities for these systems. Their 
approach, based on large effective core potentials com- 
bined with appropriately set core-polarization potentials, 
was particularly successful in predicting binding energies 
and spectroscopic properties of the alkali-metal dimers in 
the ground and low-lying excited states. 

In this paper we report a systematic ab initio study 
of the isotropic van der Waals Cg coefficients for the 
alkali-metal atom -f molecule (A-I-AB) and molecule -|- 
molecule (AB -I- AB) systems. We also derive charac- 
teristics essential for simple models of the corresponding 
ultra-low-energy collisions. In the following calculations 
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masses of the bosonic '^Li, "^^Na, ''^K, ^'^Rb, and ^^^Cs 
isotopes were used. 

II. METHODOLOGY 

The purpose of this paper is to provide essential pa- 
rameters for modelhng of collisions between the polar 
molecules in their ground rovibrational state. If collid- 
ing molecules are in j = states only isotropic part of 
the interaction potential governs its scattering properties 
at very long range larger than i?vdw = {'^mCa/h'^)^ , 
where m is the reduced mass of the colliding system. If 
the strength of the anisotropy of the potential becomes 
comparable with spacing of the appropriate rotational 
energy levels of the molecule, then anisotropic term be- 
comes important: for example in case of A-I-AB collisions, 
the coupling driven by Cq2 between j = and j = 2 
channels becomes important if the potential anisotropy 
is comparable with 6B. The same argument holds also 
for the AB+AB collisions. 

It is well known [28, 29] that within the Born- 
Oppenheimer approximation the isotropic van der Waals 
Cq coefficient of a two-monomer system (X-t-Y) can be 
decomposed into two contributions each corresponding to 
a different effect. First, it contains a dispersion contri- 
bution Cq^^^ that physically represents the interaction of 
fluctuating instantaneous dipole moments, which are due 
to the movements of electrons, which correlate between 
interacting species at long range. Secondly, in the case of 
the heteronuclear dimers, a permanent molecular dipole 
moment induces a dipole moment on the atom, which in 
turn interacts with the permanent molecular dipole mo- 
ment. This induction contribution Cg"*^ is usually smaller 
than the dispersion contribution. 

The dispersion contribution to the isotropic van der 
Waals Ce coefficient can be calculated from the following 
integral 

o /■oo 

= - / ax(iw)aY(«w)dw (1) 
Jo 

where i is the unit imaginary number, uj is frequency, and 
Q;moi(i<^) = -^["xxliw) + ayy{iuj) + a^^{iuj)] (2) 

is the orientation-averaged molecular dynamic dipole po- 

larizability. The indiiction contribution to the isotropic 
van der Waals Cf, coefficient can be expressed as 

C^^ = MiaY(0)+M^ax(0), (3) 

where is the corresponding permanent molecular dipole 
moment. If the monomer X is an atom in the spherical- 
symmetry ground state and the overlap of the charge dis- 
tribution of interacting species can be neglected, the first 
term in the above equation vanishes. The total isotropic 
Cg coefficient of the atom + molecule system is a sum of 
the dispersion and induction contributions. 



For molecule + molecule systems in their ground rota- 
tional state there also exists a (non-Born-Oppcnheimer) 
rotational contribution to the effective isotropic Cg re- 
sulting from a second-order coupling of the dipole-dipole 
term [30-32]. It which has the form C^°' = ix^/QB where 
B is the molecule rotational constant. Then the total 
isotropic Cg coefficient for the molecule -|- molecule sys- 
tem is a sum of the dispersion, induction, and rotational 
contributions. 

Proper choice of the electron basis set is crucial for 
quantum chemistry calculations of the dipole moments 
and polarizabilitics. For lithium and sodium atoms we 
have used available core-valence correlation-consistent 
basis sets cc-pCV5Z designed by Prascher et al. [33], 
which we augmented by one set of diffuse functions. 
Effective-core potentials (ECPs) with tailored valence ba- 
sis sets for heavy (K - Fr) alkali metal atoms have been 
optimized by Lim et al. [34]. These ECPs are small- 
core type potentials, i.e. the outermost 9 electrons are 
described explicitly. To eliminate possible errors due to 
the basis incompleteness we have improved the original 
valence basis sets by adding g and h functions, and aug- 
menting the basis sets by one set of diffuse functions. 
These basis sets have been tested on the atomic static 
dipolar polarizabilitics, which have been calculated with 
the spin-restricted open-shell coupled cluster method [35] 
with single, double and non-iterative triple excitations 
[RCCSD(T)] employing a finite-field approach. In all 
cases the agreement with the reference values of Dere- 
vianko et al. [36] was very good (the difference for Na 
was 2.2 a.u.; less than one atomic unit for other alkali 
metals). 

All alkali-metal dimers in their ground electronic state 
X^5]+ (near their equilibrium lengths) have their excited 
states significantly separated in energy, thus we can prop- 
erly describe them by a single-reference Slater determi- 
nant, which is ideal for using the coupled cluster ap- 
proach [35]. For the molecular calculations we took the 
equilibrium distances, which were optimized by the Paris 
group [26, 27]. 

For the molecular dynamic polarizability calcula- 
tions we employed the time-independent coupled clus- 
ter polarization propagator method in singles- and dou- 
bles approximation (TI-CCSD). This was introduced by 
Moszynski et al. [37] and implemented in Molpro 2010.2 
program [38]. Several approximations to the full time- 
independent polarization propagator were discussed by 
Korona et al. [39]. In our study, we used the so-called 
CCSD(3) approximation of the TI-CCSD method, which 
is exact to the third order of the electronic correlation op- 
erator. In benchmark calculations against the dynamic 
dipole polarizabilitics based on the full-configuration- 
interaction response functions, the CCSD(3) approxima- 
tion demonstrated systematically a smaller error than the 
other approximations introduced there [37, 39]. Finally, 
in this paper we have used a finite-field CCSD(T) ap- 
proach in order to evaluate dipole moments and static 
dipole polarizabilitics of the alkali-metal dimers. Such 
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calculations were also needed to verify the accuracy of 
the TI-CCSD dynamic polarizabilities. 



III. RESULTS AND DISCUSSION 

The dipolc moments of the heteronuclear alkali-metal 
dimers calculated with the finite- filed CCSD(T) method 
as the first derivatives of energy with respect to the elec- 
tric field applied arc collected in Table I). The FF- 
CCSD(T) data are in a good agreement with those ob- 
tained by Aymar and Dulieu [26], with an error of at most 
10% for LiNa and KRb molecules - note that for these 
species the dipole moment and charge separation between 
atoms is significantly smaller than in other cases and the 
corresponding dipole moments are small. The dipole mo- 
ment we have obtained for KRb (0.62 D) is somewhat 
larger than the experimental value (0.566±0.017 D [4]) 
and the value in Ref. [32]. For the sake of consistency 
we used our calculated values in further calculations. In 
addition to the finite-field CCSD(T) values we also calcu- 
lated the finite-field CCSD values in order to check how 
important is the inclusion of triply excited configurations 
in calculations of the alkali-metal dimer dipole moments. 
The FF-CCSD values are also in good agreement with 
our reference FF-CCSD (T) data (the FF-CCSD values 
are systematically higher by 10%). The fact that the 
triples contribution to the dipole moments is not too sub- 
stantial indicates that the FF-CCSD(T) result might be 
very close to real values as the expansion of the molecu- 
lar wavefunction in terms of number of excitations should 
converge rather quickly. 

In order to verify the quality of the molecular dynamic 
dipole polarizabilities calculated with TI-CCSD we per- 
formed further tests by checking their values in the static 
limit against the polarizabilities calculated with the FF- 
CCSD(T) approach and literature data. As the refer- 
ence values we have used those published by Deiglmayr 
et al. [27] who used a 2-electron full configuration in- 
teraction method with carefully tailored large-core effec- 
tive core potentials including core polarization potentials. 
This approach has proven to be accurate, for example, in 
predicting experimental values of the dipole moments of 
KRb [4] , LiCs [40] , and transition dipole moments RbCs 
[41] . We have also used for comparison the values of Ur- 
ban and Sadlej [42], which were obtained with entirely 
different approach - using all-electron basis sets with 
scalar relativistic effects included. Our finite-field results 
agree very well with the results from Refs. [27, 42]; our 
FF-CCSD(T) values of the orientation-averaged molecu- 
lar static dipole polarizability are systematically right in 
between their values with differences not exceeding 6%. 
The agreement between our TI-CCSD and FF-CCSD (T) 
values of the orientation-averaged molecular static dipole 
polarizability is even better (see Table I). With the ex- 
ception of LiNa, where the difference is indeed negligible, 
the TI-CCSD values are systematically higher than the 
FF-CCSD (T) values with the differences never exceeding 



2.5%. The anisotropy of the molecular static dipole po- 
larizability Aamoi(O) exhibits the same tendency. With 
the exception of LiNa the TI-CCSD values are system- 
atically higher than the FF-CCSD(T) values with the 
differences ranging from 2.5% to 13%. 

In the evaluation of formulas (1) and (3) we used the 
TI-CCSD values of the molecular dynamic dipole polariz- 
abilities, the FF-CCSD(T) values of the molecular dipole 
moments and molecular static dipole polarizabilities, and 
the values of the atomic static and dynamic polarizabili- 
ties from Ref. [36]. The integral in Eq. (1) was evaluated 
using Gauss quadrature for 50 frequencies also provided 
by Derevianko et al. [36]. 

Tables II and III contain the predicted isotropic van 
der Waals Ce coefficients for the A+AB and AB-I-AB 
systems, respectively. In the case of the A-I-AB systems 
there is a very clear progression in increase of the Cq 
coefficient for both A and B from Li toward Cs. The 
induction contribution to Cq is usually small; only in the 
case of significantly polar LiCs, LiRb, NaCs, and NaRb 
molecules it is within the range 10-23%. Our Cq values 
for the K-hKRb, Rb-hKRb, Rb-hRbCs, and Cs-hRbCs 
systems are systematically larger than those reported by 
Kotochigova [24] by 8%, 6%, 35%, and 41%.. respectively. 
The result for KRb-|-atom is clearly in good agreement 
with the result of Kotochigova, however, the difference 
for RbCs-fatoms is significantly larger. The dynamic 
polarizability in Ref. [24] has been obtained as a sum- 
over-statc with appropriate transition dipole moments of 
the RbCs molecule. It is likely that this way the RbCs 
dynamic polarizability might have been underestimated 
using such procedure by neglecting some contributions or 
underestimating the continuum contribution. Note also, 
that the induction contribution for atom-hdiatom has not 
been included in Ref. [24]. 

It was proposed recently to approximate the Cq co- 
efficients by simply adding the pairwise atom-atom Cq 
coefficients [24, 43]. Our calculations have verified this 
model as seemingly reasonably good for heavy atoms (Cs 
and Rb) interacting with weakly polar molecules. This 
nice agreement is, however, fortuitous, since this additive 
approximation includes only dispersion and no induction 
. In this case, the approximation of the trimer dispersion 
forces by simply adding them among dimers overcompen- 
sates the lack of the induction interaction. 

As expected, the effective isotropic Cq coefficients for 
the AB-I-AB systems very strongly depend on the AB 
dipole moment. Only the LiNa and KRb dimers with 
smallest dipole moments are dominated by electronic 
contribution to the Cq coefficient, in other cases rota- 
tional contribution completely dominates the total Cq 
coefficient. For the AB-I-AB systems, there is also a very 
distinct pattern in increase of the electronic contribution 
similar to the A+AB systems. The KRb isotropic Cq 
coefficient is higher by 10% compared to the value of Ko- 
tochigova [24] and by 6% with respect to the value given 
by Buchachcnko and coworkers [32]. Our results are in 
agreement with those reported by Quemener et al. [25] 
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TABLE I. Dipole moments (in Debye), rotational constant (in cm~^), the orientation-averaged molecular static dipole polariz- 
ability (in atomic units), and the anisotropy of the molecular static dipole polarizability (in atomic units) of the ground X^E"*" 
states of hctcronuclcar alkali-metal dimcrs calculated at the equilibrium interatomic distances from Ref. [26]. 



dimer 


re/ao [26] 


B/hc 




amol(O) 


amoi(O) 






LiNa 


5.4518 


0.425 


0.48 


237.7 


237.6 


156.3 


155.7 


LiK 


6.268 


0.293 


3.41 


324.2 


326.9 


234.5 


240.7 


LiRb 


6.5 


0.254 


3.99 


347.2 


352.1 


262.0 


272.7 


LiCs 


6.93 


0.218 


5.39 


391.9 


399.1 


317.8 


333.1 


NaK 


6.61 


0.094 


2.72 


358.1 


362.7 


247.2 


260.9 


NaRb 


6.88 


0.070 


3.31 


387.1 


393.9 


279.2 


299.7 


NaCs 


7.27 


0.058 


4.63 


439.3 


448.0 


339.4 


364.1 


KRb 


7.688 


0.037 


0.62 


523.5 


532.3 


367.6 


409.5 


KCs 


8.095 


0.030 


1.98 


596.0 


606.8 


436.1 


488.9 


RbCs 


8.366 


0.017 


1.32 


638.6 


653.0 


462.1 


531.1 



TABLE II. The isotropic Cs van der Waals coefficients (in 
atomic units) for the alkali-metal A -|- AB systems. The last 
column shows the value based on pair-wise atom-atom addi- 
tive model. 



atom 


dimer 


y^disp 




Ce 




Li 


LiNa 
LiK 
LiRb 
LiCs 


2217 
2885 
3098 
3452 


6 

294 
407 
740 


2223 
3179 
3505 
4192 


2856 
3711 
3934 
4454 


Na 


LiNa 
NaK 
NaRb 
NaCs 


2358 
3405 
3673 
4092 


6 
187 
275 
539 


2364 
3592 
3948 
4631 


3023 
4003 
4239 
4783 


K 


LiK 

NaK 

KRb 

KCs 


4821 
5364 

7428 
8298 


520 
334 
17 
175 


5341 
5698 

7445 
8473 


6219 
6344 
8171 
9056 


Rb 


LiRb 

NaRb 

KRb 

RbCs 


5688 
6357 
8154 

9751 


790 
539 
19 

87 


6478 
6896 
8173 
9838 


7235 
7373 
8964 
10353 


Cs 


LiCs 
NaCs 
KCs 
RbCs 


7652 
8555 
10995 
11772 


1803 
1324 
242 
110 


9455 
9879 
11237 
11882 


9911 
10073 
12005 
12509 



for the LiNa (difference of 4%), LiK (20%), LiRb (17%), 
and LiCs (11%) systems. Thcisc values are very sensi- 
tive to the dipole moment and rotational constant of the 
molecule, thus even small differences in these character- 
istics can easily translate into a 20% difference in the 
dominating rotational part of the Cg coefficient. 

Known Cg coefficients allow us to determine the energy 
limits for single partial-wave scattering. The p-wave or 
d-wave scattering starts to dominate if the collision en- 
ergy is comparable to the appropriate centrifugal barrier 
heights: for the A-I-AB collisions it is the p-wave, while 
for the bosonic AB-f AB collisions it is the d-wave scatter- 
ing. In Table IV we have included the centrifugal barrier 
heights for the; A+AB scattering and bosonic AB-I-AB 
scattering. Their values approximately determine the 
single partial-wave regime. The same table contains also 
the mean scattering lengths [44], which illustrate, in a 
sense, a characteristic length scale of the corresponding 
interaction potential. 



IV. CONCLUSIONS 



TABLE III. The isotropic Ce van der Waals coefficients (in 
atomic units) for the alkali-metal AB-I-AB systems. 



dimer 


y^disp 




CI"' 




LiNa 

LiK 

LiRb 

LiCs 


3582 
6024 
6963 

8670 


17 
1167 
1711 

3520 


110 
404491 
876031 

3397216 


3709 
411682 
884705 

3409406 


NaK 

NaRb 

NaCs 


7461 
8696 
10822 


820 
1313 
2916 


508325 
1497080 
6932958 


516606 
1507089 
6946696 


KRb 

KCs 


14202 
17716 


62 

723 


3456 

450681 


17720 

469120 


RbCs 


20301 


345 


160336 


180982 



In conclusion, we have reported a complete ah initio 
study of the isotropic Cg van der Waals coefficients in 
all possible A-I-AB and AB+AB systems, where A and B 
are two distinct alkali-metal atoms and AB are molecules 
are in their ground state. Given the rapid development of 
the field and many ongoing experiments with polar alkali- 
metal molecules, we expect that these results will be ben- 
eficial for modelling their coUisional properties, which 
are crucial for stability studies of the ultracold molec- 
ular dipolar gases in traps. In future studies we would 
like to pay increased attention to the role of anisotropy 
in ultracold collisions, and we would also like to exploit 
our results when constructing potential energy surfaces 
for various coUisional systems. 
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TABLE IV. Heights (in /x/f) of p-wave centrifugal barriers 
Vp for the A+AB systems and d-wave centrifugal barriers Va 
for the AB+AB systems, respectively, with mean scattering 

leii,y,tlis a (iu <ui) for the eorrespoudiiis collisions. 
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atom 


dimer 


Vp 


a 


dimer 


K-i 


a 


Li 


LiNa 


2442 


39 


LiNa 


2293 


30 




LiK 


1844 


44 


LiK 


111 


110 




LiRb 


1581 


46 


LiRb 


28 


157 




LiCs 


1397 


48 


LiCs 


8 


244 


Na 


LiNa 


684 


49 










NaK 


380 


58 


NaK 


64 


125 




NaRb 


300 


61 


NaRb 


16 


187 




NaCs 


256 


65 


NaCs 


6 


251 


K 


LiK 


221 


68 










NaK 


177 


71 










KKd 


liz 




KKd 


iiy 


04 




KCs 


95 


85 


KCs 


15 


155 


Rb 


LiRb 


64 


86 










NaRb 


56 


89 










KRb 


47 


95 










RbCs 


32 


104 


RbCs 


17 


131 


Cs 


LiCs 


29 


105 










NaCs 


26 


108 










KCs 


23 


112 










RbCs 


19 


117 
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